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ABSTRACT 


This report represents an annual progress report describing, 
rather briefly, the research conducted on this grant. In those cases 
where the findings are already available in the open literature, 
the material is merely highlighted and appropriate references given. 

There are five phases to the project which are described sepa- 
rately. The general subject matter of the phases includes: 

1. Ignition of fuel drops by a shock wave and passage of a shock 
wave over a burning drop. 

2. The energy release pattern of a two-phase detonation with 
controlled drop sizes. 

3. The attenuation of shock and detonation waves passing over 
an acoustic liner. 

4. Experimental and theoretical studies of film detonations. 

5. A simplified analytical model of a rotating two-phase 
detonation wave in a rocket motor 
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I. INTRODUCTION 


The research covered by this third annual progress report 
represents a continuation of our efforts devoted to the study of de- 
tonation waves in liquid-gas systems. The motivation for the work 
is associated with liquid propellant rocket motor combustion in- 
stability although certainly the studies are also applicable to in- 
ternal combustion engines, jet propulsion engines, safety aspects 
of spilled liquid fuel, coal mine explosions, and weaponry. The re- 
search has been divided into 5 phases, although all of them are in- 
timately related. For the most part these phases are briefly summar- 
ized and the reader is referred to other publications for a more 
complete treatment. The exception to this is where the material 
herein represents the only printed information available on the 
particular facet of the problem. 

Phase A has been primarily concerned with the breakup and 
ignition of fuel drops by shock waves. The experimental portion 
of this study as well as a theoretical treatment of the ignition 
behavior was completed in the past year. The research is now 
concentrating on the passage of a shock wave over a burning drop. 

Phase B has been devoted to the assessment of the approximate 
energy release pattern in two phase detonations insofar as they affect 
the significant overpressures observed. 



Emphasis has been placed on detecting the variations observed 
when 2 discrete drop sizes, as compared to a monodisperse system, 
are utilized. Attention has also been given to an analytical prediction 
of the overpressures which, along with the experimental results, 
should yield the major energy release pattern. 

Phase C has been aimed at the feasibility of arresting two phase 
detonations or the development of shock waves into a detonation by 
means of an acoustic liner. The engine case with a rotating detona- 
tion wave as well as a shock tube configuration have been considered. 

Phase D has been concerned with film detonations, that is where 
the liquid fuel is in the form of a thin film on the wall of a tube. The 
experimental work and the two dimensional analysis were completed. 
The effort on this phase is now limited to preparing a NASA report 
on the subject. 

Phase E consists of an analytical study of a rotating two phase 
detonation wave in an annular rocket motor. This problem has been 
essentially completed although a modest effort to include the radial 
effects is underway. 

The progress made on each of the 5 phases will now be described. 
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II. R ESEAR CH R ES ULTS 


Phase A - Shock Wave Ignition of Fuel Drops 
Predictions of Wake Characteristics at Ignition. 

Earlier in the program a number of experiments were performed 
wherein shock waves of varying strength were passed over liquid fuel 
drops in oxidizing as well as inert gaseous atmospheres. The aero- 
dynamic shattering and ignition characteristics for diethylcyclohexane 
and n-hexadecane drops in oxygen and oxygen- nitrogen mixtures were 
determined for varying incident Mach numbers, different ambient pres- 
sures, and for differing initial drop diameters. Immediately after col- 
lision of the shock with the drop, liquid fuel, in the form of very small 
drops, is carried into the turbulent wake of the parent drop. Here the 
droplets vaporize and mix with the surrounding gases. After an ig- 
nition delay time, provided conditions are sufficient, ignition will occur 
in the wake. Experimentally it has been found that the ignition delay 
time depends on the fuel, initial drop diameter, the Mach number of 
the incident shock, and the initial oxidizer partial pressure. The 
significance of the Mach number dependence is primarily associated 
with the fact that the shock Mach number determines the dynamic 
pressure and temperature of the convective flow over the drop. The 
results of these experimental investigations are reported in detail 
in references 1, 2, and 3. 
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In order to gather some information on the conditions required in 
the wake of the drop before ignition could occur, an analysis was con- 
ducted which utilized some experimental information. This analysis 
is described in detail by Kauffman 3 . For the purpose at hand the 
analysis will just be outlined and a few representative results shown. 
First of all the analysis assumes a single drop and utilizes experi- 
mental results to describe the trajectory of this parent drop. The 
convective flow behind the shock and relative to the drop is then used 
in a boundary layer stripping model to calculate the size and rate of 
addition of microspray droplets to the wake region. The predicted 
variation of miscrospray drop size with time as the parent drop 
flattens transverse to the gas flow and breaks up is shown in Fig. 1. 

As indicated on the figure, the results are shown for a 1520 p di- 
ethylcyclohexane drop exposed to a Mach 3. 75 shock in 20 in. Hg 
of pure oxygen. The predicted microspray drop sizes are small. 
Further, the analysis predicts that this size varies with the square 
root of the parent drop diameter. It is realized that the boundary 
layer stripping model used is not necessarily the dominant mechanism 
throughout the entire breakup history, particularly in the later stages. 
However, it was considered adequate for the purpose at hand and was 
used for lack of anything better. The variation of the peak micromist 
size for a range of shock strengths and initial pressures is shown in 
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Fig. 2. The decrease in size with Mg and/or with increased initial 
pressure is primarily due to the importance of dynamic pressure. 

Combining the foregoing with the experimental observation that 
the microspray drops are almost instantaneously accelerated to the 
local gas velocity, the rate of evaporation of these drops at the local 
static pressure and temperature was calculated (admitting that the 
appropriate equation is not known for such dynamic conditions and 
under the influence of surface tension). As a result of all of the mech- 
anisms involved , the composition of the wake will not be uniform in 
distance downstream of the parent drop. A calculation of such a 
variation is shown in Fig. 3 wherein the vaporized fuel- oxidizer ratio 
is shown as a function of distance behind the forward stagnation point 
of the parent drop. The variations for different times are shown for 
a 1520jU DECH drop upon interaction with an M = 3. 50 shock wave in 20 
in, Hg of 02- The experimentally observed ignition time was 77 ps. 
The stoichoimetric f/o ratio is 0. 292 by mass so that the maximum 
f/o ratio achieved in the wake is approximately 1/3 stoichiometric. 

The lack of detailed knowledge of the kinetics of DECH oxidation 
has precluded the attainment of an ignition criteria derived from first 
principles. Accordingly, it has been necessary to examine the cal- 
culated state of the wake at the experimentally determined time of 
ignition; as was just done in Fig. 3. 

The calculated total mass of fuel in the vapor phase in the wake 
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at ignition for a 15 20 ju drop subjected to various strength shocks at 
3 different initial pressures is shown in Fig. 4. The plot indicates 
that for the reduced pressures of oxygen appreciably richer f/o 
ratios are required for ignition at any of the Mach numbers cal- 
culated and for reduced Mach numbers at a fixed pressure. The 
latter characteristic is presumably a temperature effect. However, 
the one atmosphere case showed essentially no variation with M s . 
Similar plots made for other drop sizes indicate that more fuel 
vapor is in the wake at the onset of ignition for the larger drop 
sizes. This is borne out by the experimental observation that 
stronger blast waves are experienced with the larger fuel drops. 

Finally, it is of interest to examine the sensitivity of the results 
to the uncertainty in the micromist size. Fig. 5 shows the variation 
of the maximum value of f/o in the wake at ignition for various shock 
Mach numbers and for the microspray drop size taken at 1/2, 1. 0, 
and 2. 0 times the liquid boundary layer thickness. The range ex- 
tends over an equivalence ratio of about 1/8 - 2. While this spread 
is appreciable, it is encouraging that the predictions are at least 
reasonable. 

Interaction of a Shock Wave with a Burning Drop. 

The foregoing has been concerned with the ignition of a fuel 
drop by a shock wave. However, in actual combustion devices fuel 
drops are often burning before interaction with a pressure pulse. 
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Accordingly, an experimental study of the interaction of a shock 
wave with an initially burning fuel drop has been started. 

The primary objective of this phase of the study is to determine 
the alteration of burning rate of the drop when exposed to a step 
pressure increase. 

The apparatus which injects and ignites a single fuel drop into 
the shock tube is described in this section with some of the prelimi- 
nary results of shock-drop interaction experiment. 

The method of ignited drop injection is based on the sudden retrac 
tion of the assembly of a hypodermic needle, from which the drop 
is initially suspended, and an electrode, which gives an electric 
spark, shortly before ignition. The drop injection system is shown 
in Fig. 6. It essentially consists of a needle- electrode assembly 
and spring loaded plunger which is locked by the locking pin prior 
to the run. The solenoid retracts the pin and the plunger pulls up 
the needle- electrode assembly by a thin steel wire. The timing for 
activating the solenoid and spark electrode is accomplished by the 
shock wave through a pressure switch which is located 8. 052 ft. 
upstream of the drop injection point. The ignition spark goes be- 
fore the solenoid unlocks the spring loaded plunger, although the 
same signal from the pressure switch is used. The drop injection 
system is covered by a dome to facilitate evacuation and filling of 
the shock tube. 
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There are two major difficulties with this drop injection system 
which are not yet completely solved. First, the attainment of a 
single spherical droplet is difficult and secondly, there are certain 
technical limits in shortening the time interval between shock wave 
sensing and complete retraction of the needle- electrode assembly 
out of the shock tube. 

The needle is made from a 0. 035 in. o. d. x 0. 022 in. i. d. tube 
with a 0. 020 in. o. d. tube fitted inside. 

Fig. 7 shows drops obtained from this system (non ignited, 
numbers indicate the time elapsed from the instant of short cir- 
cuiting the pressure switch). During many trials to get accept- 
able shape of the drops, the following observations were made: 

(a) The quantity of the fuel suspended before retraction should 
be a maximum. If the quantity of the fuel suspended is too 
small, the drop shape is distorted. 

(b) The lateral movement of the needle during retraction should 
be minimized; otherwise, the drop is thrown laterally, dis- 
torted, and sometimes disintegrated into small droplets. 

Fig. 7 shows a little lateral movement of the needle in 
spite of considerable efforts to eliminate it. 

(c) The suspended drop is distorted by the surface tension force 
when the needle is suddenly retracted. This makes the 
shape of the injected drop like a tadpole. 
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In hopes of eliminating the tail of the "tadpole”, steps are 
being taken to minimize the opening of the fuel flow area at the 
needle tip. The use of a smaller diameter needle is also planned 
for obtaining more spherically shaped drops. 

The time interval allowed between the instant of switching by 
the shock wave and the instant of shock wave-drop interaction is 
extremely small. For instance, at shock Mach number M s = 2. 0, 
the time interval is 3750 / us and at M g = 3. 0, 2500 jus. It is found 
that the major sources of prolonging retraction of the needle- elec- 
trode assembly are the current rise time of the solenoid and associ- 
ated movement of the locking pin and the acceleration of the mass 
of the plunger and the needle- electrode assembly. The solenoid, 
which is rated at 24 VDC, is operated by a high voltage, short 
duration current discharge from condensers. Earlier, the con- 
densers were charged to 320 V and the unlocking time was about 
1800 ps. This time decreased to 1300 ps when the voltage was in- 
creased to 580 V. The unlocking time is determined from a dis- 
continuity in the current-time trace from the oscilloscope. 

The sum of the mass of needle- electrode assembly and plunger 
is about 15 grams. The spring constant of the spring used is 0. 49 kg/cm 
and the initial upward force to the plunger is 8. 4 kg. Simple analysis 
of a spring and mass system is in good agreement with the experi- 
ment. Experiment and the analysis shows that it takes about 3500 ps 


9 


for the needle to travel 10 mm upward. 

Since there are no obvious ways of modifying the present electro- 
mechanical system to appreciably shorten the time interval, efforts 
are being directed to activate the appropriate time delay units from 
the time of bursting of the shock tube diaphragm. 

At the higher shock Mach numbers, the needle- electrode assembly 
is still in the shock tube at the time of shock arrival and is damaged. 
Some effort is being made to devise a suitable protection of the needle 
electrode assembly without causing unfavorable disturbance to the 
interaction phenomena. 

A series of preliminary experiments using fairly weak shock has 
been made. Fig. 8 shows the sequence of pictures of interaction of 
an ignited drop with shock waves. Although the disturbance from 
droplet shape, i. e. from the tale of the ’’tadpole", is considerable, 
most of the features of interest of the interaction are retained. The 
shock Mach number is 1. 55 and the nominal drop diameter is 2040 ji- 
lt is noted from Fig. 8 that the shock is accelerated drastically in 
the hot region. The initial laminar flame is blown off the drop by 
the convective flow behind the shock. The blown off flame apparent- 
ly ignites the micromist in the wake region in between 50 to 70 ps 
after initial interaction. Fig. 9 shows the sequence of the interaction 
of a shock wave with a drop which has not been ignited. The break up 
characteristics are not noticeably different from those of the ignited 
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drop. However, it is seen from the comparison of Fig. 8 with Fig. 9, 
that the micromist in the wake region is consumed much more in the 
case of the initially ignited drop. 

In these cases where the flame attaches to the wake, the subse- 
quent shed micromist will be continously consumed and hence the 
severe blast wave phenomenon, experienced in shock wave ignition^ 
will be precluded. 

Accordingly, it is worth while to determine the condition under 
which the flame will remain attached to the drop, attached to the 
wake, or completely blown off. To this end, streak schlieren 
photography will be used in experiments utilizing various shock 
strengths and initial oxygen concentrations. The affect on burn- 
ing rate will also be established. 


11 



Phase B - Energy Release Patterns 


Four series of experimental tests were conducted subsequent to 
the previous report^. Two of these involved polydisperse spray deto- 
nations consisting of DECH as fuel in a pure gaseous oxygen environ- 
ment, while the other two were concerned with the dilution of the charge 
gas with unreactive and reactive gaseous additives. 

As was done previously, the composition of a polydisperse spray 
was described in terms of the total equivalence ratio, and the 
partial equivalence ratios, of the constituent drops, where _ 

23 Also, the "equivalence ratio fraction" of species i is a^ = 
i 

<#>i/ 0T- 

In the first series of tests, sprays were detonated which consist- 
ed of 300 p and 1300 p diameter drops, with a^QQ = 0. 525 and = 0. 80. 
This bimodal distribution represents the maximum separation in drop 
sizes possible using the polydisperse drop generator^, under the 
condition of approximately equal division in total available energy 
between the two species. A monodisperse spray of 1300 p diameter 
drops was detonated for comparison at = 0. 80, but the corres- 
ponding monodisperse spray of 300p diameter drops could not be 
produced at this equivalence ratio, since 38 drop-producing cap- 
pilaries would have been required of the drop generator which has 
a 20-needle capacity. 


12 



The structural differences between the monodisperse (1300p) spray 
and the bimodal (1300/i, 300 jll) spray detonations can be observed in 
Fig. 10. The apparent ignition plane of the 1300/n drops (at about 2 in. 
behind the leading shock) is more distinct in the monodisperse than 
in the bimodal case. The intensity of the explosive ignition of the large 
drops has been somewhat subdued by the earlier ignition of the small 
drops. In addition, the mean pressure level is more sustained in the 
monodisperse than in the bimodal spray. This is because, in the bimodal 
spray, the small drops, once ignited, add energy in the region ahead of 
the large drop ignition plane, causing the pressure level there to drop; 
in the monodisperse spray of large drops, the pressure remains fairly 
constant until the drops ignite. 

The average detonation velocity in a monodisperse spray of 1300/i 

drops at <f> =0.8 was U = 5220 ft/sec while in the bimodal spray it was 

X s 

U = 5890 ft/sec, an increase of about 11%. This results in a lower aver- 
s 

age pressure in the monodisperse case (even though it is more sustained) 
by about 25%. It appears that with a large drop size difference in a poly- 
disperse spray, the factor which most greatly influences the average pres- 
sure levels in the wave is the extent to which shifts in distribution affect 
wave velocity and the leading shock wave pressure. With lesser drop size 
differences, the wave velocity is not as greatly affected by distribution, and 
in that case, the average pressure levels appear to be related to distribu- 
tion through the tendency towards increasingly sustained pressures as the 
proportion of large drops is increased. 
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Another series of experiments examined the effect of total equiv- 
alence ratio on detonation properties. Bimodal sprays of 300pand 750 ju 
diameter drops were used. The tests were conducted at constant equiv- 
alence-ratio fraction, 0. 550, with total equivalence ratios vary- 

ing between = 0. 265 and <p^ = 0. 775. 

The decrease in droplet spacing accompanying an increasing 0,p was 
expected to result in progressively smoother energy release. That is, 
due to the increasingly probable proximity of small drops to the wakes 
of the large drops, the explosive ignition of the former could be expected 
to cause ignition in these wakes earlier than normal. This should re- 
sult in a decrease in the intensity of large-drop blast waves. 

The expected tendency of decreasing overpressure ratio with in- 
creasing was observed (Fig. 11). However, the pressure structure 
smoothness did not improve over the range of </>rj. studied. 

Although the peak overpressures in Fig. 11 are decreasing, the 
magnitude of the average reaction zone pressure increases with in- 
creasing $ (Fig. 12). This is due to increasing shock pressures 
resulting from an increasing total heat addition, which is also accom- 
panied by an increasing wave velocity (Fig. 13). 

The third series of tests dealt with the effect of charge gas dilu- 
tion by nitrogen. Monodispersed DECH sprays consisting of 750 ju 
diameter drops were employed in these experiments, with the charge 
gas mixtures at .1 atm pressure. 
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The equivalence ratio of DECH with respect to the oxygen present 
in the charge gas was held at <£rp = 0. 914 in all tests. 

The detonation velocity as a function of nitrogen mole fraction, 
X N is plotted in Fig. 14. Nitrogen mole fraction is the ratio of 
the number of moles of nitrogen to the total number of moles of 
nitrogen and oxygen in the charge gas. The greatest^ content 
which produced detonation was X-^ = 0. 5. Prior to this point the 
detonation velocity appears to be only slightly affected by the 


presence of the N 2 - 

The average pressure level in these waves (Fig. 15) decreases 
significantly with increasing nitrogen concentration. This is due to 
a decreased intensity of the explosive ignitions by individual drops 
in the spray, which reduces the average overpressures. Peak overpressures 
were not as strongly affected by N 2 content. 

When the diluting additive is a reactive gas, the detonation 
system is of mixed type, consisting of gas-phase and two-phase 
parts. These are referred to as "hybrid" detonations. In such detonation, 
the gas-phase reaction is presumably completed well before the drops 
ignite. Hence the equivalence ratio of the gas-phase component is 
defined in terms of the initial oxygen present. On the other hand, the 
equivalence ratio of the two-phase part may be expressed in terms of 
either the initial oxygen concentration, or the oxygen concentration 
remaining after completion of the gas-phase part. In the experiments 
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presently described, the equivalence ratio of the liquid fuel is 
measured with respect to the initial oxygen content. These tests 
were conducted in monodisperse sprays of 480p diameter drops, with 
hydrogen gas as the diluent. The equivalence ratio of DECH was held 
at unity. 

Wave velocity and Mach number as functions of hydrogen con- 
centration are shown in Fig. 16, along with the corresponding gas- 
phase hydrogen-oxygen curves for comparison. At low hydrogen 
mole fractions, the two-phase detonation supports the H 2 -O 2 reaction, 
whereas at high hydrogen concentrations, the H 2 -O 2 detonation supports 
the two-phase part. 

Very pronounced structural changes occur with increasing hy- 
drogen content. These are shown on Fig. 17. The apparent ignition 
zone position of the liquid drops moves rearward in a continuous 
manner. Moreover, the addition of heat to the flow field in the 
region just behind the leading shock by the H2-0 2 reaction produces 
progressively reduced overall reaction zone levels. 

After X H = 0. 5, droplet ignition is so delayed that there exist in 

Lt 

the records from these waves extended regions of relatively minor 
pressure excursions behind the leading shock and the voilently react- 
ed liquid drops appear to act as a ’’piston" in preventing the trailing 
rarefaction behind the gas detonation. After Xj^ = • 667 (stoichio- 
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metric), no overpressures appear, and the fuel from the liquid drops 
probably reacts much more gently. 

Many of the resutls just described were presented in a paper 
entitled, "Two-Phase Detonations with Bimodal Drop Distributions, " 
which was presented at Third International Colloquium on Gasdynamics 
of Explosions and Reactive Systems, Marseille, France. The paper 
has been accepted for publication in Astronautica Acta. 

The theoretical description of the reaction zone structure in a 
monodisperse two-phase spray detonation has progressed. The igni- 
tion time estimates described previously^ have been compared with 
revised experimental data (Fig. 18 and 19) showing improved correla- 
tion, particularly in the case of the "whole wake" thermal imbalance 
theory. This theory predicts the occurrence of explosive ignition 
when the rate of heat production due to pre- ignition reaction in the 
wake region of a given parent drop in the spray exceeds the rate of 
energy removal from this region due to turbulent mixing with the ex- 
ternal stream. 

A detailed analytical account of the time- varying processes which 
take place in the reaction zone of a two-phase detonation has been pre- 
pared in this analysis, each individual droplet in a fuel spray through 
which the detonation passes is presumed to undergo the sequence of 
events described by the "whole wake" ignition theory. The blast waves 
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produced by the individual wake explosions colalesce, after a certain 
distance, into two planar shocks, one of which moves towards the 
leading shock (forward-moving) and the other which moves away from 
it (rearward- moving). 

The model is constructed in a one-dimensional coordinate system 
which moves at the average detonation velocity. Fuel droplets arrang- 
ed in a square array pass through the leading shock, and after travel- 
ing a distance into the reaction zone corresponding to the overall igni- 
tion delay time, a given plane of drops ignites. The details of the three- 
dimensional process by means of which the individual spherical blast 
waves (from each drop in this plane) coalesce into planar waves is 
not analyzed. However, the time from ignition to coalescence, dur- 
ing which three-dimensional effects are important, is assumed to be 
negligibly small. The one-dimensional distributions of the gasdynamic 
variables within the blast wave region just after coalescence occurs is 
forced to satisfy integral conservation of mass, momentum, and energy 
with respect to that which existed in the same physical region just 
prior to ignition. 

Thereafter, the one-dimensional blast wave is allowed to freely 
expand within the reaction zone. For this purpose, a two-step (pre- 
dictor-corrector) numerical scheme^ is utilized to solve the Navier- 
Stokes equations within the blast wave region. 

Eventually, the forward- moving shock from the blast wave will 


18 



overtake the leading shock, imparting to it an increased velocity. This 
is the mechanism of wave sustenance; i. e. , between successive colli- 
sions of this type, the leading shock decays. Hence the detonation does 
not actually propagate at a constant velocity, but oscillates about an 
average speed. 

The detonation structure contains a timewise periodicity resulting 
from the regularly occurring explosive ignitions. The frequency of 
these ignitions is coupled to the average wave velocity and the spacing 
between the droplet planes. After each ignition, the newly formed 
blast wave itself expands within the expansion wave of its predecessor. 
Several such nested blast waves can be present within the reaction 
zone at a given instant. 

Only a fraction of the total energy available in a drop is releas- 
ed when it ignites. A significant mass of fuel survives in the form 
of microspray and in the parent drop. Subsequent to ignition the drop- 
shattering and microspray evaporation processes continue, with the 
difference that it now occurs in a much hotter and more turbulent 
region. 

For purposes of analysis, the post- ignition reaction process 
is decoupled from its complex aero thermo dynamic interaction 
with the convective flow, and a simplified energy release rate is 
assumed. It is also assumed that the drops remain in their origi- 
nal orientations so that each plane of drops acts as a separate, 
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moving heat source within the flow in the reaction zone, downstream 
of the ignition point. 

The results from this calculation are graphically displayed 
utilizing a high speed (1200 baud information exchange rate) Computek 
CRT terminal. Plots of pressure, density, temperature, and gas 
velocity are recorded as functions of position within the wave, at 
given times in the calculation. In addition, a computer-animated 
film follows the time- varying structure of the detonation through 
several cycles. 

An example of the data generated by this calculation is present- 
ed in Fig. 20. The abscissa of the plot denotes distance from the 
explosion center; i. e. , the position within the reaction zone at 
which the periodic explosive ignitions occur. In producing the 
movie, this plot would represent one image, of which four frames 
on the film would be taken. 

When this calculation is completed on a given set of conditions, 
there is effectively available a field of data in the x-t plane (within 
limits on x and t) within which all information is known, including 
the shock positions and strengths and the values of the gasdynamic 
variables between shocks. Therefore, additional information can 
be gained by simulating the movement of a pressure transducer 
through this field. 
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In leading shock fixed coordinates, such a transducer pene- 
trates the reaction zone along a path in the x-t plane described 
by dx/dt = U„ (the detonation velocity). At any given instant, the 

o 

output from the transducer represents an integration of the pressure 
distribution over its exposed surface. 

For a circular transducer of 0. 1 in. diameter, a typical trans- 
ducer "sweep" of this type (through the data from which Fig 20 was 
taken) is shown on Fig. 21. The abscissa is now transducer distance 
behind the leading shock. 

It can be seen that the shocks appear to the transducer as round- 
ed "bumps". There is also a rise time even at the leading shock, 
corresponding to the time of leading shock passage over the trans- 
ducer. 

The one dimensional-time unsteady-structure of a two phase 
detonation just described will be presented in a paper for the forth- 
coming 14th International Symposium on Combustion, The Combustion 
Institute, Pennsylvania State University, August 1972. 

A detailed coverage of all of the material of this Phase B will 
appear shortly in the Ph. D. thesis of T. H. Pierce. 
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Phase C - Acoustic Liner Studies 


The study of the effectiveness of acoustic liners as attenuation 
devices for two -phase detonation waves has been the principal aim for 
this phase of our work. Fully recognizing that to be effective the acoustic 
liner has to damp out compression waves generated in the combustion 
zone and thereby hopefully preclude the development of the combustion 
(or flame) front into a detonation, different aspects of the interaction of 
shock waves and detonation waves with acoustic liners have been studied. 

So far the limited experimental observations and the approximate analyt- 
ical calculations have indicated that for relatively weak shocks of strength 
comparable to the pressure pulses in the combustion zone the shock 

velocity reduction obtainable by interaction with an acoustic liner is 
6 *7 

relatively small ’ . The mass efflux from behind the shock wave into 
the liner cavities appears to be the main cause of the velocity change 
in the shock wave. 

An analysis to study the effect of mass bleed behind the leading shock 
of a two phase detonation wave has been completed. Jump conditions for 
a two phase mixture of burned gases and un burned fuel and oxidizer 
droplets with different incoming velocities were developed which employed 
the following main assumptions (Fig. 22): 

1. An adiabatic and frictionless wall, 

2. A one -dimensional C-J two -phase detonation propagating at 
constant velocity. 

3. The mass efflux and the momentum loss are expressed as fractions 
of the corresponding values of the incoming gases and droplets. 
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An F-number analysis of the type first used by Adamson and Morrison 
was carried out to yield the relations for the properties of a C-J detona- 
tion with mass bleed. 

By including additional conditions of cyclic behavior and by linking 

the average pressure on the injection plate to the design equilibrium 
4 

pressure > the above analysis was extended to cover the case of an 
annular chamber with a rotating two -phase detonation wave with mass 
bleed. The results of calculations in typical cases of a tube detonation 
and a rocket motor wave are shown in Fig. 23 and 24. A parametric 
study of the two phase wave with mass bleed, in the two cases of an 
annular chamber and the laboratory tube conditions, indicates the ways 
of correlating the experimental results in the laboratory to engine condi- 
tions without actually resorting to motor testing. Numerical results 
indicate that in order to obtain any appreciable velocity reduction for the 
wave the gas bleed fraction has to be large. 

Estimates of the mass efflux behind the leading shock of the detonation 
were made using different wall porosity values. Drop breakup dynamics 
were included in the estimation of reaction zone length. Thus, by combin- 
ing these estimates with the parametric calculations of effect of bleed, 
the steady state velocity of the wave with bleed can be determined. After 
considering numerous calculations of the above type using typical liner 
open-area ratio, and drop size values, the indications are that for a 
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steady propagating wave in the motor case it is not possible to obtain signifi- 
cant reductions in propagation velocity of the wave. This is mainly due 
to two factors: (i) the available reaction zone length is too small to obtain 
any significant mass bleed and (ii) the 'effective open -area' ratio of the 
liner is small. The situation is somewhat better in the case of typical 
laboratory test conditions compared to that in the engine chamber. 

After reaching the above conclusions the next step taken was to study 
the effect of gas mass bleed on the blast waves originating from the drop 
combustion zone and moving towards the leading shock front. If these 
blast waves are sufficiently weakened by the gas bleed before they catch 
up with the leading shock, the sustaining mechanism for the wave propaga- 
tion is disturbed, thereby causing the leading shock to slow down. A 
one -dimensional time -dependent finite difference technique was used 
with drop size, leading shock strength and wall porosity as parameters. 
Early results indicate that for a well developed detonation the blast wave 
strength is not appreciably reduced in the time it moves through the reac- 
tion zone. Correlations between the engine chamber and laboratory test 
conditions are to be established. 

Based on the analytical work done to date, conditions for the labora- 
tory tests on the effects of bleed on two-phase detonations, correspond- 
ing to chamber wave parameters when the liner effectiveness is optimum, 
will be chosen and the attenuation effectiveness of the acoustic liner will 
be experimentally observed. 


24 



Phase D - Film Detonations 

During the past year Dr. C. S. R. Rao completed the thesis 
"Theoretical and Experimental Study of Film Detonations” and was 
awarded the Ph. D. Degree. Key parts of this thesis appear in the 
paper "A Two Dimensional Theory for Two Phase Detonation of 
Liquid Films, " by C. S. R. Rao, M. Sicnel, and J. A. Nicholls 
which was prepared for publication and appears in Combustion 
Science and Technology , Vol. 4, pp 209-220, 1972. In this paper 
and Dr. Rao's thesis the simple one dimensional film detonation 
theory, which is reported in "A Simple Theory for the Propagation 
of Film Detonations, ” by M. Sichel, C. S. R. Rao, and J. A. Nicholls, 
Thir teenth Symposium (International) On Combustion , pp 1141-1149, 
(1970), has been extended to take into account the effect of boundary 
displacement on the flow in the reaction zone. 

The analysis is based on the model shown in Fig. 25 below, which 
shows the reaction zone of a typical film detonation from coordinates 
moving with the wave. In these coordinates the shock is stationary 
and the wall or surface of the shock tube moves with the propagation 
velocity U . A turbulent boundary layer forms immediately behind 

D 

the leading shock and results in drag forces and heat transfer from 
the core of the reaction zone to the liquid film causing the liquid 
to vaporize. A key hypothesis in the analysis is that vaporization 
sets the rate at which the fuel in the film reacts with the oxidizer 
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in the core of the reaction zone. Support for this hypothesis comes 
from the agreement between theory and experiment as indicated 
below. The main difference between the one dimensional and two 
dimensional theories is that in the latter the displacement effect of 
the boundary layer upon the core flow is taken into account. 

The propagation characteristics of the film detonations are 
obtained from the mass, momentum and energy conservation equa- 
tions across the detonation front with the effects of drag, heat 
transfer, and boundary layer displacement within the reaction zone 
included in the formulation. The end of the reaction zone on the 
Chapman Jouquet (C-J) plane is taken as the point at which the 
film is completely vaporized. A detailed analysis of the Chapman- 
Jouguet condition shows that this choice for the C-J plane is valid 
as long as the fuel air ratio is not much larger than the stoichio- 
metric value. Combustion is assumed to occur within a diffusion 
flame sheet imbedded in the boundary layer. The boundary layer 
parameters are determined by combining the results of shock tube 
boundary layer theory with relations which have been developed in 
the analysis of combustion in hybrid (liquid-solid fueled) rockets. 

Theoretically computed propagation speeds, reaction zone 
lengths, and pressure ratios are compared to experimentally mea- 
sured values in Figures 26, 27, 28 which summarize the key results 
of this investigation. The theory and measurements shown in 
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Figures 26, 27, 28 are for film detonations propagating down a 
square tube with internal dimensions of 1. 64 in. x 1. 64 in. filled 
with oxygen at atmospheric pressure and temperature and with di- 
ethylcyclohexane as the fuel. Film detonation with one, two and 
four walls of the shock tube wetted with fuel were considered. In 
these figures L is the reaction zone lengths, P 3 /P 1 is the pressure 

ratio across the detonation and u s /u so is the ratio of the actual 
propagation velocity u g to the propagation velocity u so of a pre- 
mixed gaseous detonation without losses and with the same fuel and 
oxidizer. Results from both the one and two dimensional theories 
are shown. 

Figures 26, 27, 28 indicate good agreement between the analysis 
and experimental measurements and support the model of the film 
detonation described above. In particular it appears that vaporiza- 
tion is the rate limiting process, at least in the case considered here 
with very thin films and overall fuel air ratios of the order of the 
stoichiometric value. 

The effect of diluting the oxygen with nitrogen has also been 
investigated. Theory and experiment are compared in Figures 
29, 30, 31 for oxygen mass fractions of 1. 0, . 89, and 0. 675. 

Again agreement between theory and experiment was good. In the 
experiments it was found that film detonations would not propagate 
when the dilution exceeded a certain critical value. Thus it was im- 


27 



possible to detonate diethylcyclohexane films in air. The conclusions 
may, of course, change for different fuels, particularly fuels with a 
high vapor pressure. Also, higher initiation energies may lead to 
steady detonation in air. 

At the present time a contractor report describing the details of 
the film detonation analysis is under preparation. The study of film 
detonation will then be terminated. 
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Phase E - Theoretical Analysis of a Tangential Two- Phase De tona tion 

This phase of the research has been devoted to a theoretical treat- 
ment of a tangential mode strong discontinuity type of instability. The 
major portion of the work was completed and written up in the past year. 
A detailed treatment of the work appears as the Ph. D. thesis of I-wu 
Shen, 1971, entitled, "Theoretical Analysis of a Rotating Two-Phase 
Detonation in a Liquid Propellant Rocket Motor. " A paper, based on 
this work with the same title, was presented at the Third Internation- 
al Colloquium on Gas Dynamics of Explosions and Reactive Systems, 
University of Marseille, Marseille, France, September, 1971. This 
paper will appear soon in Astronautica Acta. A short summary of the 
work completed and that in progress will now be given. 

The analysis assumes that a one-dimensional two-phase detona- 
tion propagates at some constant, but unknown, angular velocity in a 
thin annular combustion chamber. Thus, no radial variations are 
taken into account. The walls are taken to be adiabatic and friction- 
less. It is assumed that the drops remain at their initial condition 
until they are converted to gases in the thin reaction zone. The 
detonation is treated as a discontinuity and the flow between waves 
is assumed isentropic with no interaction between drops and the burn- 
ed gases. Further, the influence of the wave on propellent injection 
rate and combustion efficiency is neglected. The rotating wave is 
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assumed to decay to a sound wave as the nozzle is approached. 

In this problem, as posed, there are 9 variables inasmuch 
as the wave velocity is unknown a priori. The variables in addition 
to wave velocity are the pressure, density, temperature and gas 
velocity immediately upstream as well as immediately downstream 
of the wave. In order to determine these 9 unknown, 9 equations are 
needed. These are obtained by employing the jump conditions for a 
two-phase mixture of burned gases, unburned fuel and oxidizer drop- 
lets. The cyclic condition provides additional information and, finally, 
the pressure distribution along the injector plate (determined by the 
method of characteristics) is integrated numerically to obtain the 
average pressure which is then linked to the design equilibrium 
chamber pressure. The analysis then predicts the variation of wave 
strength (as measured by pressure ratio or wave velocity and assum- 
ing the wave exists) with nozzle area ratio, specific heat ratio, in- 
jection velocity, impingement distance, equilibrium chamber sound 
speed, chamber diameter, and fuel distribution. 

It is found that the pressure ratio across the wave may be reduced 
(i. e. stabilizing) by reducing the contraction ratio, reducing the pro- 
pellant flux density at the outer periphery near the injection plate, 
increasing the injection velocity, reducing the chamber diameter, and 
decreasing the chamber speed of sound. The wave strength is indepen- 
dent of chamber length and propellant mass flow. An increase in the 


30 



latter increases the pressure level of the engine but the pressure 
ratio across the wave remains the same. 

A further portion of the analysis has considered the influence of 
drop size insofar as it affects the validity of the analysis, and hence 
the possibility of realizing such an instability. 

A comparison of the predicted results with the few available 
experimental results are very encouraging, both as to the pressure 
predictions and limits of validity of the analysis . 9 

A NASA report streamlining the above analysis for the designer's 
purpose is being prepared. Also, a small effort is continuing to in- 
clude the first effects of radial variations such as would exist in a 
cylindrical chamber. The equations have been set-up for this addition 
but predicted results are not yet available. 
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Diameter of Micromist vs Dimensionless Time, 1520 /u DECH Drops. 
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Fig. 2. Maximum Diameter of Micromist vs Shock Wave 
Mach Number, 1520M DECH Drops. 
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Fig. 5. Maximum Value of Fuel to Oxidizer Ratio vs Incident 
Shock Wave Mach Number, Various Micromist Sizes. 
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Figure 6. Ignited Drop Injection System. 
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Fig. 7. Drop Formation; Non-Burning. 
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Fig. 8. Shock Wave Interaction of Burning DECH Drops. M = 1. 55, d = 2040 /z. 
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Fig. 10. Pressure Distribution behind Detonation; $rp = . 8 , 
Pi = 1 atm; Upper: Monodisperse Spray, 1 300 ju Diameter 
Drops; Lower: Bimodal Spray, 300 ju and 1300 p Diameter 

Drops, a t 300 = . 5 25. 
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Fig. 12. Average Pressure During Initial 30 psec Period After Passage of 
Leading Shock in Spray of 300juand 750/i Diameter Drops vs Equivalence 

Ratio; a 75 o = 0. 55, Pj = 1 atm. 





Fig. 13. Detonation Velocity vs Equivalence-Ratio in Spray 
of 300juand 750p Diameter Drops; a^gQ = . 55. 
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Fig. 14. Detonation Velocity vs Nitrogen Mole Fraction; 
D 0 = 750] u, Pi = 1 atm, <f> = 0. 914. 
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Fig. 18. Comparison of Theoretical Ignition Time with Experimental Data, 
Whole-Wake Thermal Theory; DECH Drops in Pure 02, D 0 = 930jU (Experi- 
mental Data from Ref. 3). 
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Notation for 
Fig. 22-24: 


( ) - gas, ( ),- drops, (1) - behind wave, (s) - shock, 
0 (o) - ahead of wave 


f = bleed fraction 

f = fraction of V 0 the bleed mass 
m ® 

has in leaving C. V. 

= droplet mass flux fraction 
Md = (Pd v d)/(Po V o + ^d v d) 

Q d * = (h d - h Q 0 )/(CpT 0 ) = heat released 
by droplets 

Q * = <h Q 0 " h l 0) / <C P 0 To) = ^ 

s released by gases 
m = mean molecular weight 
y = ratio of specific heats 


ER = equivalence ratio 

K d = V V o 

p l*= p l/ p o 
p l* = p l /p o 
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M = V / d 
o o' a 

P * = P /P 
c c' o 

T * = T /T 

C C O 

V * = V /d 
w w' c 


Fig. 22. Jump Conditions for C-J Detonation with Mass Bleed. 
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Fig. 25. Two-Dimensional Model for Film Detonations 
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Fig. 26. Propagation Characteristics of One Wall Film Detonations — 
Comparison of Experiment and Theories 
(Fuel - Diethylcyclohexane; Oxidizer - Oxygen) 
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Fig. 27. Propagation Characteristics of Two Wall Film Detonations 
Comparison of Experiment and Theories 
(Fuel - Diethylcyclohexane; Oxidizer - Oxygen) 
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Fig. 28. Propagation Characteristics of Four Wall Film Detonations 
Comparison of Experiment and Two-Dimensional Theory 
(Fuel - Diethylcyclohexane; Oxidizer - Oxygen) 
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Fig. 29. Effect of Oxidizer Dilution on One Wall Film Detonations — 
Comparison of Experiment and Two-Dimensional Theory 
(Fuel - Diethylcyclohexane) 
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Eig. 30. Effect of Oxidizer Dilution on Two Wall Film Detonations — 
Comparison of Experiment and Theory 
(Fuel - Diethylcyclohexane) 
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